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Abstract. The dicyemid mesozoa, obligate symbionts in 
the cephalopod kidney, are simply organized multicellular 
animals. They have long been the subject of phylogenetic 
debates. Some authors have suggested that dicyemids rep¬ 
resent an offshoot from an early metazoan ancestor. Other 
workers considered them to be degenerated progeny of 
higher metazoa, possibly parasitic trematodes. We deter¬ 
mined the almost complete nucleotide sequences of 18S 
rDNA in two species of dicyemid, Dicvema orientate and 
Dicyema acutieephalum, isolated purely from cephalopod 
urine. We compared these sequences with sequences de¬ 
termined in the present study from three flatworm species, 
as well as with a variety of eukaryote sequences obtained 
from databases. The phylogenetic trees reconstructed with 
the use of the neighbor-joining, maximum-parsimony, 
and maximum-likelihood methods indicated that the di¬ 
cyemids belong among the triploblastic animals (Bila- 
teria). However, we cannot firmly establish the position 
of the dicyemids within the Bilateria because we cannot 
ignore the problem of long branch attraction between the 
myxozoans, dicyemids, nematodes, and acoel flatworms. 
The present results favor the hypothesis that the dicyemids 
do not represent an early divergent metazoan group, but 
rather a group degenerated from a triploblastic ancestor. 

Introduction 

The dicyemids are simply organized multicellular an¬ 
imals consisting of an outer layer of 20-40 ciliated somatic 
cells and an inner core of one long axial cell. They are 
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obligate symbionts in the kidney of cephalopods. Their 
life cycle is complex; the asexually produced vermiform 
embryos increase the population in the host, while the 
infusoriform embryos arising from fertilized eggs pass out 
of the host body with the urine and are thought to infect 
another cephalopod host (Brusca and Brusca, 1990). 

The dicyemid mesozoans have long been the subject 
of a phylogenetic controversy (Brusca and Brusca, 1990; 
Willmer, 1990). They were at first considered to be an 
extant link between the Protozoa and the Metazoa (Hy¬ 
man, 1959). However, the resemblance of their complex 
life cycles to those of parasitic trematodes has led some 
authors to propose that the dicyemids are descended from 
an established metazoan group and that their simple body 
organization results from degeneration attributable to 
parasitism (Nouvel, 1948; McConnaughey, 1951; Stun- 
kard, 1954: Ginctsinskaya, 1988). Others still view the 
simple body construction of dicyemids as truly primitive 
and hold that the group represents an offshoot from early 
divergent metazoa (Dodson, 1956; Hyman, 1959; Lapan 
and Morowitz, 1974). 

The phylogenetic relationships of eukary otes have re¬ 
cently come under intense scrutiny in the light of new 
molecular data. Phylogenetic analyses using nucleotide 
sequences of 5S rRNA suggested that a dicyemid (Di¬ 
cyema misakiense) diverged early among such lower me¬ 
tazoa as sponges, cnidarians, and flatworms (Ohama et 
at., 1984: Hori and Osawa, 1987). But phylogenetic trees 
based upon comparisons of about 120 sites in the nu¬ 
cleotide sequences of 5S rRNA were different from those 
inferred from longer nucleotide sequences of 18S or 28S 
ribosomal RNA (Field el at., 1988; Christen et ai, 1991; 
Katayama et ai, 1993; Wainright et ai, 1993; Kobayashi 
et ai, 1993). We have sequenced 18S ribosomal RNA 
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genes (18S rDNA) in two species of the dicyemid mesozoa. 
Our comparison of the nucleotide sequences of small- 
subunit rDNA for a variety of organisms indicates that 
the dicyemids belong among the triploblastic animals. 

Materials and Methods 

Biological materials 

We determined almost the entire sequence of 18S 
rDNA in two species of dicyemid mesozoans and three 
species of turbellarians (Platyhelminthes). Pure samples 
of the dicyemids Dicyemci acuticephalum and Dicvema 
orientate were collected from the urine of Octopus vulgaris 
and Sepioteuthis lessoniana, respectively (Furuya et ai , 
1992a). Specimens of Convoluta naikaiensis (Acoela) and 
Planocera multitentaculata (Polycladida) were collected 
on the shore near the Ushimado Marine Laboratory. 
Specimens of Dugesia japonica (Trieladida) were obtained 
from the brook near the Ushimado Marine Laboratory. 
All were frozen quickly and kept at — 80°C until use. 

In addition to the sequences of the above five species, 
we used the sequences of 23 eukaryotes—including ani¬ 
mals, protists, plants, and fungi—for which almost com¬ 
plete 18S rDNA sequences were available in databases. 
The species used and their accession numbers are as fol¬ 
lows: Paramecium tetraurelia (Ciliophora), X03772; 
Oxytricha nova (Ciliophora), X03948; Crypthecodinium 
cohnii (Dinozoa), M64245; Theileria annulata (Apicom- 
plexa), M64243; Sareocystis muris (Apicomplexa), 
M64244; Uartmanella vermiformis (Rhizopoda), 
M95168; Saceharoinyces cerevisiae (Fungi), JO 1353; Fi- 
lobasidiclla neoformans (Fungi), X60183; Arabidopsis 
thaliana (Plantae), X16077; Volvox carteri (Plantae), 
X53904; Berne eueumis (Ctenophora), D15068; Tricho- 
plax adhaerens (Placozoa), LI0828; Anemonia sulcata 
(Cnidaria), X53498; Scvpha ciliata (Porifera), L,10827; 
Jlenneguva sp. (Myxozoa), U13826; Caenorhabditis ele- 
gans (Nematoda), X03680; Afoliniformis moliniformis 
(Acanthocephala), Z19562; Schistosoma mansoni (Tre- 
matoda), X53047; Crassostrea gigas (Bivalvia), X60315; 
Artemia salina (Crustacea), X01723; Sagitta crassa 
(Chaetognatha), D14363; Asterias amurensis (Asteroidea), 
D14358; and Xenopus laevis (Vertebrata). X04025. 

DNA isolation 

Genomic DNA was extracted by the method described 
previously (Wada et ai, 1992). In brief, the frozen samples 
were lysed in TE buffer (10 m M Tris-HCl, 0.1 AI EDTA, 
pH 8.0) containing 0.5% sodium dodccyl sulfate. After 
digestion with proteinase K (100 pg/nA) at 50°C for 3 h, 
DNA was extracted with phenol and precipitated in 
ethanol and an equal volume of 5.0 M ammonium ace¬ 
tate. Samples resuspended in TE buffer were further pu¬ 


rified by RNase A digestion (20 pg/m\) at 37°C for 1 h 
followed by ethanol precipitation. 

Amplification of 1SS rDNA 

The 18S rDNA was amplified by the polymerase chain 
reaction (PCR; Saiki et ai, 1988) in an Air Thermo-cycler 
1645 (Idaho Technology). Almost the entire length of 18S 
rDNA was amplified using synthetic oligonucleotides, 
5 -CTGGTTGATCCTGCCAG-3' (primer 0) and 5'- 
CCTTGTTACGACTT-3' (primer 10) as the terminal 
primers. Amplifications were performed in 50 ^1 of 
50 m M Tris-HCl (pH 8.5), 250 pg/m\ BSA, 2 m M Mg 2+ , 
with 0.2 m l/ each dNTP, 50 p M primers, template DNA 
(5-10 ng), and 2U Taq DNA polymerase (TOYOBO). 
The temperature regimen for 35 cycles was 20 s at 94°C, 
30 s at 50°C, and 90 s at 74°C. 

Determination of DNA sequences 

After purification of the amplified DNA by electro¬ 
phoresis in a 0.8% agarose gel, the nucleotide sequence 
was directly determined by dideoxy chain-termination 
(Sanger et ai, 1977) using Sequenase ver 2.0 (USB) and 
[ 35 S]-dATP (Amersham). All DNA samples were se¬ 
quenced in both directions and from several separate am¬ 
plifications with terminal primers (0 and 10) and internal 
primers. The internal primers used were primer-1 (5- 
CCGGAGAGGGAGCCTGA-3'), primer-2 (antisense of 
primer-1), primer-3 (5-CAGCAGCCGCGGTAATT-3'), 
primer-4 (antisense of primer-3), primer-5 (5'-GCGAA- 
AGCATTTGCCAA-3'), primer-6 (antisense of primer 5). 
primer-7 (5'-GAAACT(TC)AAAGGAAT-3'), primer-8 
(antisense of primer-7), and primer-9 (5-ACGGGC- 
GGTGTGT(AG)C-3'). The positions corresponding to 
these primers in 18S rDNA sequences are shown in Figure 
1. The continuity of the DNA fragments was confirmed 
by overlapping of the sequences. 

Phylogenetic analyses 

Sequences were aligned manually on the basis of max¬ 
imum nucleotide similarity (Fig. 1). Alignment gaps were 
inserted to account for putative length differences between 
sequences. Some regions could be confidently aligned and 
were presumed to be homologous. However, we could 
not unequivocally determine the optimal alignment for 
the regions containing deletions, insertions, or highly 
variable sequences. We excluded positions from the anal¬ 
ysis according to the following rule: Positions where a gap 
was present for any taxon were not used in analyses. In 
the present study, we found this simple rule alone ade¬ 
quate for excluding the regions of ambiguous homology 
(the regions where two or more equally optimal align¬ 
ments were present) from the analysis because in those 
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regions alignment gaps were always serially inserted in 
many sequences. The phylogenetic trees were recon¬ 
structed using the PH YL1P package version 3.5c (Fclsen- 
stcin, 1989) and fastDNAML (Olsen et al 1993). Tree¬ 
building procedures used were the neighbor-joining (Sai- 
tou and Nek 1987), the maximum-parsimony (Fitch, 
1971), and the maximum-likelihood (Felsenstein, 1981). 
For the neighbor-joining analysis, evolutionary distance 
values were calculated by the formula of Jukes and Cantor 
(1969). The degree of support for internal branches of the 
trees in the neighbor-joining and the maximum-parsi¬ 
mony trees was assessed by bootstrap levels of support 
(Felsenstein, 1985) determined by 500 bootstrap repeti¬ 
tions. 

Results 

In D. acuticcphahim, D. orientate, Convoluta naikaien- 
sis, Dugesia japonica, and Planocera multitentacnlata, al¬ 
most the entire length of 18S rDNA was amplified by 
PCR from the genomic DNA. The sequence (1500— 
1700 bp) was determined directly from PCR products. 
The sequences have been deposited in databases (GSDB, 
DDBJ, EMBL, and NCBI) under the following accession 
numbers: D26529 for D. acuticcphahim; D26530, D. ori- 
cntalc; D17558, Convoluta naikaiensis; D17560, Dugesia 
japonica; D17562, Planoccra multitcntaculata. To infer 
the phylogenetic position of the dicyemids within the eu¬ 
karyotes, we aligned the almost complete nucleotide se¬ 
quences of 18S rDNA of the above five species with the 
23 eukaryote sequences we obtained from databases. Me¬ 
tazoan taxa were chosen to represent phyla broadly; pro¬ 
tozoan taxa were chosen to represent the more recently 
derived groups. Figure 1 shows a sample of the alignment 
for 9 out of 28 species included in the present analysis. 
This alignment reveals that throughout the eukaryotes 
the sequences are highly conserved in some regions and 
highly variable in others. After exclusion of the regions 
of ambiguous homology, 1070 sites (Fig. 1) remained for 
phylogenetic inference. Phylogenetic trees shown were re¬ 
constructed by the neighbor-joining (Fig. 2). the maxi¬ 
mum-parsimony (Fig. 3), and the maximum-likelihood 
(Fig. 4) analyses. 

Among the phylogenetic trees reconstructed by the 
three methods, the topologies were largely congruent with 
one another, though branching with low bootstrap support 
within the metazoan lineage showed somewhat conflicting 
arrangements. The metazoans—including triploblasts 
(Bilatcria), diploblasts, dicyemids, and a mvxozoan— 
formed a monophylctic assemblage in the three trees. 
Within the metazoan assemblage, triploblasts formed a 
discrete monophyletic unit together with the mesozoa and 
the Mvxozoa. The branches of triploblasts were in general 
longer than those of other taxa. 


The grouping of the dicyemids with the triploblastic 
animals was supported by a bootstrap value of 100% in 
both the neighbor-joining (Fig. 2) and the maximum-par- 
simonv (Fig. 3) analyses (because of the enormous com¬ 
putation time required, bootstrapping was not performed 
in the maximum-likelihood analysis). In the trees recon¬ 
structed by the three methods, dicyemids were grouped 
with Cacnorhahditis elegans (a nematode), Henneguya 
sp. (a mvxozoan), and Convoluta naikaiensis (an acoel 
flatworm), though bootstrap confidence level for this 
grouping was low. 

To corroborate the inclusion of the dicyemids in the 
triploblastic lineage, we analyzed subsets of taxa shown 
in the present paper as well as several different sets of taxa 
including some of the following species (the accession 
numbers for 18S rDNA data are shown in parentheses): 
Cryptomonas phi (X57162), Babesia bovis (M87566), 
Tctilla japonica (D15067), Svcon calcaravis (D15066), 
Mncmiopsis leiclyi (L10826), Tripcdalia cystophora 
(L10829), Paraspadclla gotoi (D14362), Antedon serrata 
(D14357), Strongylocentrotus intermedins (D14365), 
Balanoglossus carnosus (D14359), Oikoplcura sp. 
(D14360), Branchiostoma jloridae (M19571), Homo sa¬ 
piens (X03205). In all sets of taxa analyzed, triploblasts 
formed a monophyletic unit and the dicyemids were 
placed within the triploblastic clade with high bootstrap 
confidence level. Grouping of Caenorhabditis, Convoluta, 
Henneguya, and Dicyema were consistently observed. 
Trichoplax, another enigmatic animal whose phylogenetic 
position is controversial (Brusca and Brusca, 1990; Will- 
mer, 1990), was always positioned outside the triploblastic 
assemblage, confirming analyses by Wainwright et al. 
(1993). 

Discussion 

The present molecular phylogenetic study based upon 
comparisons of nucleotide sequences of 18S rDNA shows 
that triploblastic animals form a monophyletic assemblage 
within the metazoan subtree and that the dicyemid me¬ 
sozoa are an ingroup of the monophylctic unit of triplo¬ 
blastic animals. Monophyly of triploblastic animals has 
repeatedly been shown in previous molecular phylogenetic 
analyses with 18S or 28S rDNA sequences (Field et al., 
1988; Christen et al. , 1991; Wainright et al., 1993; Ko- 
bayashi et al. . 1993; Smother et al., 1994). The present 
analysis confirms the recent claim by Smother et al. (1994) 
based upon 18S rDNA sequences that the Mvxozoa are 
closely related to the triploblastic animals. 

As in phylogenetic trees previously constructed on the 
basis of rDNA sequences (Christen et al., 1991; Wainright 
et al., 1993; Smother et al., 1994), Trichoplax, which had 
once tentatively been grouped in the phylum Mesozoa 
(see Brusca and Brusca, 1990), was positioned within the 
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D.acuticephalum ATGCATGCGTAAGCTCATGCTCT-CTAATGAGCGAAACCGCAGACGGCTCATTAAATCGGACATAACTTACTTGATAATCTC-ACACGACGGTGGATAAC 
D. orientale nnnnnn.TA. .. .T.-.G.-... 
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D.acuticephalum 
D.orientale 

C. naikaiensis 

D. japonica 
P.multitentaculata 
A.sulcata 

P.tetraurelia 
A .thaliana 
S.cerevisiae 


TGTTGTAAATCTAGAGCTAATACATGCGTACAGCTTCTCG-AAGCGCAATTATTAGA- 
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C. .G_T.A . G . . TCC . GACTTCT- -GGAAGGGAT .T.T.TTCAAAACCAATGCGGGTTCT- 

C. .G_T.C . A .TACCGGACGC AAGAATGGTA. . .T.TTTAA--CCAT- 
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C. .G_T.T. . A.ATC . . GACCCTTT-GGAAG .GAT. T.T.TAAAAAATCAATG-T—CTTCG- 
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A.sulcata 
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-ATTAGTTGACTCTGGATAAACTTGCTGATCACA-GnAGTGCCTGCGACGTATTAATT-AGATGTCCGCCCCATCAATTT 
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Primer 1,(2) 
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Figure 1. Nucleotide base sequences of 18S rDNA from two dicyemids, Dicycma acuticephalum and 
Dicyema orientale. and three flatworms, Convoluta naikaiensis. Dugesia japonica, and Planocera multiten- 
taculata The sequences are aligned to that of small-subunit rDNA of a diploblasl. Anemonia sulcata; a 
prolist, Paramecium tetraurelia; a plant. Arabidopsis thaliana. and a fungus. Sacchromyces cerevisiae. A 
period indicates that the base at that position is identical to that in D. acuticephalum. a hyphen indicates a 
gap. and an “n" indicates an undetermined site. The 1070 positions used for phylogenetic inference are 
shown by single lines above the alignment. The positions corresponding to the internal primers are shown 
by double lines below the alignment. The primers whose numbers are in parentheses are antisense to the 
sequence shown here. 
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.T.A.A.CT.T. .GCCCG—. .TG. CCGGTC. . A. TT . TTCGT . T-AC . . GAT . TCCAACG . GGC . . T . CCTTCTGGCTA 
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750 800 

-GC--TTTACCTTGAACAAA-ATAGAGTGCTTAAGGC-AAGCATTCTGCTTG-A 

-GCTACG. .-...-A.A. 

CATCGATGATATTAGTTGCACTTTGTTGTGACGACTAGTGGAAACGGTGTT..C..T....GT...-T.T.C..A..-..C.T..GG.A...C. 

GTTAAGTGCACTTTATTGGGATCTTTTACAA--TAACCGACAA-G.T.A. . .-T.A . »G . G . . T . ATGCT . GCAT 

TTTTGGTGCTCTTAATTGAGTGCCTTAATTGCCCG-GCCA-CG.T.A. . .AT.G.C. .AT.-.G. .CCAA. . .C. .A. 

CGCGTGTGCTCTTGACTGAGTGTGCGCGGGAGTTGCGA-CG.T.A. . .-T.T.C. .A. .-.G. .-CAGC.A. 

-TTAGGGTTGCAGCTGGGCGAGTAG-ACAA..A. . .-T.T.CC--.G.T-C.TC. .CG.A. 

CGCTCCTGGTCTTAATTGGCCGGGT-CGTGCCTCCGGCG-CTG. . . .T.G . . .-T.C. .A. . .C.A.GCTC. .G. 

ACCTTGAGTCCTTG—TGGCTCT-TG-GCGAACCAGGA-CT.T.A. . .-T.T .C . . A . . - . G. . G . AT-C .A. 


_ _ 850 _ __ 900 

D.acuticephalum TATCTAAGCATGGAATAATAGAATAAGAC-TTTTCTA-TTGGTT-ACG-A-TAGTAAAAGTAATG-TTAACAGAGACAGCCGGGGGCATCCGT 

D .orient ale .C.A.-.G--T..C-.G.A.A.T... 


C.naikaiensis 

. .ATAA. 

...G..GA.T. 

..-C.GAG.CAACGTTTTGAACGCGGT.GCTGTTGCTCG. 
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AT GTT 
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.A. 

. . .A. 
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C.GA.CTCG..CC... 
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D.acuticephalum ATTGCTCCGTTATAGGTGAAATTCGTAGATCGGTGCAGGACGTA-CTACAGCGAA-GCATTTGCCAAC-ATGTTTTCATTAATCAAGAACGACAGTTGGA 
D.orientale .A...C..G.A.-.-.G-. 


C. naikaiensis 

D. japonica 
P.multitentaculata , 

A.sulcata 

P.tetraurelia 

A.thaliana 

S.cerevisiae 








ggtggt c g 

. .TG. . 
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.T .G. 

...ATC...A. 
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..GA. 
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. _T .G . 
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..TTA..A.A. 
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Primer 5,(6) 

_____ 1050 _ _ 1100 

D.acuticephalurn GTATCGAAGTCGATTAGATCTCGTCGTAGTTCCAACTATAAACGATGCCAACTAGCACTCCGCTGATAATATTTGGATACATTAGCGGGTGGCTTCCGGG 
D.orientale .A-G..--TTTG.T. 

C. naikaiensis _TA. . .AG.AC.C.A. . .T. . .C.TT.A.TA.GA.TT.TCCC .GTCCC- .CTCG . . TGGGCAA . AATTTAA. 

D. japonica .G.A....C....AC_C.TG. .C.T.T. . GA.GAGTTAGC . . AGGATA . .TA. .CATC .T. .TA.AA.TCA. 

P.mul titentaculata .GT.A... .C .... AC . . ..c.TG. .C.G...A....T..CG.T. GCGA. TTCGATCC . A. . . .CA. .C. 

A. sulcata .GC.A_C_AC_C.AT. . .C.G.GGA. .A.AGAG.GT. . ..GACC.CTTT. . CAC . . .AG. . . 

P. tetraurelia .G. . .A. . .A_C_AC.CTT.T. .A. .G.AG. . .G . AA .GGT . ATA- . ATTAGTCCCTTTC . . CAT .G . AA . A . 

A. thaliana .GC.A_C_AC_C_CT-C.G. .C. .GGA. . A . . G . . . GT . GC . .AT.GGAC.CC. . T . . CAC . . . AT . A . 

S. cere vis i ae . G...... . AT ...C....AC.CTT ...C.T.G.... —GA. . G. G . . G . GT . T. . .TA . . GACCC . CTC ...AC...A..A. 
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_ _ 1150 _ 1200 

D.acuticephalum AAACCAAAGTCTTTGGGTTCCGGGGGAAGT-ATGGTTGCAAAGCTGAAACTTAAAGGAATTGAC-GGAAGGGCACCACCAGGAGTGGAGCTTGCG-CTCA 
D. orient ale .-.-.-.-.-.-....C. 

C. naikaiensis .T.A.T.-.ATC.-.C.......T. 

D. japonic a .A.G.-.. .-C — — . .T. 

P.multitentaculata .-.T....A........A.C....G..T. 

A . sulcata .T.T.G...-.-.C-- .T. 

P. tetraurelia ...T.T.G.. .-.C-G.-.*.C-G. .T. 

A. thaliana . . .T.T.G. . .-.C....G.-.C-G. .T. 

S. cerevisiae . . . T.T.G. .. -.C .. ..G.. . ..C . ...G.. 

Primer 7(8) 


1250 __ 13 0 0 

D.acuticephalum ATTTGACTCAACGCAGAAAAACTCACCCGGGCCGAACACAGTGAGGATTGACAGACTGATAGCTTTTTCTTGATACTGTGGGTAGTGGTGCATGGCCGTT 
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_ 1 3 5 0 _ 1400 

D.acuticephalurn -GTTGGTTGGTTGGAGTGATTTGTCTGGTTTATTCCGATAACGAACGAGACTCAGACTTACTTACTAGTCGTACGATTCTAGTTGT-ACGAC 


D.orientale 
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TC, 

.CA. . 







.A.A.. 


GCA.TTGCTG- 

- GTTAT-.C. . 









1450 




1500 


D.acuticephalum TACTAAGAAGGATCAGTG-TGAAAACACTTGAAAATGAGCAATAACAGGTCTGTGATTGCCCTTAGA-CGTTCGGGGC-GCACGCGTGCTAC 

D. orient ale ..,-T...G---.TT. . .ACAC.T.-.-.T. A.C. 

C. naikaiensis .T. .T. . . G . C . CTGT . CCGAATTTAAAT .CGGA . AGC . . . GTGAG.-.-T....T..A.T.C. 

D. japonic a .T..T.A. .A. A. A-GCG.CT.CGT.AA.A.-.A. . . ,-T. .C.GC . ..C. 

P.multi tentaculata .T. .T. .GG . . .CA.. . .-GCAT.C. .GC. .TAC.-.-T. .C.C.C. 

A .sillcata .T..T...G.. .CTGT. .-GTG.TT. .C. . AAGTC . GGAAG.-.-T_T_C.C. 

P. tetraurelia .T. .T. . .G...CT.TGT-ATGTAAGTGCATGGAAGTTTAAG.-.C. . . .C.T.. . .C.C. 

A. thaliana .T..T...G. . .CT.TG.-CCGTT.AGGCCA.GGAAGTTTGAG.-.-T_T_C.C. 

S .cerevisiae .T. .T. . .G . . .CT.TC .-GTTTCAAGCCGATGGAAGTTTGAG.-. . ..A .... .T_.C.C. 


D.acuticephalum 
D.orientale 

C. naikaiensis 

D. japonica 

P.multi tentaculata 
A.sulcata 
P.tetraurelia 
A .thaliana 
S.cerevisiae 


_ _ _ 1550 _ 1600 

AATGAAGAAAGCAGAGGG-GGTTTTGCTTGGAAAAGCGAACTAAGCCT-TAAAATTTCTTCGTGGCAGGAATCGAGGCTTGTAATTATTCCTCG 

-G.T. . .-T. .C.C.A. G. .G.-.A. 

.C. .T.CTCTA.T.-TGAAAAACC.AAT. .A.AG.GGATTGGGA..T.ACAG..CGGAGAAATT.T.AG..G_T.CA.C.T. .G.A.. 

-GCAGTTC . .AC. A.--TAATACC. . AAC . . A. A . .T . . TTGGGC . .T.T.G.G. . . CACTG_AA . . . . G . . T .T . AA.T.A.A 

.C-T .GC .T . AC . A . --TTATTCTCC .GATCC.A. .GG.TCGGGC. .C.TGT.G. . .CC.A.A_G_G.C. . 

•C...C..TGT..AC.A.-TCTCTCC...GCC.A..GGT.TGGG...T.T.C.C...CA..G.CTG..G..A..TCA...C_T..GA..T 

.C.. .CACGTT...C.A.CTTATTTACC.G.C.CGAA.GGGCACGGGA. .T.T.G.T.GGACGTG.CTG..G. .A..TCT. ..C.AGA..T 

.C. . .T.T.TT. .AC. A.-TTCACACC .TGCC . AC . GGC .CGGG . . .T.T.-.G A.ATG. .G. .A. .TCA. . .C_G. .GG. .T 

.C. . .C.G.GC. . . C . A .-TCT AAC C . . . GCC . AG . GGT . TTGG . . .T.T.G.G. ..C.C.G.CTG. .G. .A. . . CA.G. . .T 
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A. sulcata 
P.tetraurelia 
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S.cerevisiae 


_ 1650 _ 1684 

TGCACAAGGAATTCCTAGTAATCGC-AGGTCATTAGCCTGCAATGATTACGTCCCTGCCCTTTGTACACACCGCCCGTCGCTAC 
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CA. .G.G. . .-.A.C. . .T. . .GT. . .G 


Primer(9} 
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Arabidopsis thaliana 
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Firobasidiella neofortnans 


— D icy etna acuticephalum 
Dicyetna orientale 


Diploblasts 
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Figure 2. Neighbor-joining irees showing the phylogenetic position of the dicyemids among 28 repre¬ 
sentative eukaryotic taxa. The tree was reconstructed on the basis of the pairwise distances of Jukes and 
Cantor (1969) using DNADIST and NEIGHBOR programs (PH YL1P package, version 3.5c). The tree was 
rooted by using Saccharomyces cerevisiae as an outgroup. Branch lengths are proportional to the scale given 
in substitutions per sequence position. The percentage of 500 neighbor-joining bootstrap replicates is shown 
at the node the value is supporting. 


diploblastic assemblage. This study does not support a 
close relationship between Trichoplax and mesozoans. 

These sequence data do not firmly establish the position 
of the dicyemids within the triploblastic assemblage. In 
the phylogenetic trees obtained by three different methods, 
the dicyemids formed a monophyletic unit with the 
myxozoans, nematodes, and acoel flatworms. They are 
all considered to be early divergent groups in one widely 
accepted phylogenv. The early divergence of acoel flat- 
worms in triploblastic evolution has been suggested by 
Katayama et al . (1993) from comparisons of partial 18S 
rDNA sequences. However, the myxozoans, dicyemids, 
nematodes, and acoel flatworms were all represented by 
a long branch in the phylogenetic trees (a high nucleotide 
substitution rate). Hence we cannot ignore the possibility 
that these long branches produce artifactual groupings 
within the triploblastic assemblage (Van de Peer el a /., 
1993). 

With regard to the topology of triploblast phyla, the 
present phylogenetic trees contradict those of some pre¬ 
vious analyses of 18S rDNA sequences in some points. 


The Platyhelminthes did not form a monophyletic unit 
as previously shown by Katayama el al. (1993). Mono- 
phyly of the deuterostomes and the protostomes has re¬ 
peatedly been shown in the molecular phylogenv of 18S 
rDNA (Wada and Satoh, 1994; Raff el al , 1994; Halanych 
el al, 1995), but neither group was monophyletic in all 
trees of the present analysis. Triploblastic phyla are poorly 
resolved in the molecular phylogenetic trees of 18S rDNA; 
i.e., nodes defining phyla are not supported by high 
bootstrap values. Therefore, the topology of the trees de¬ 
pends largely on the choice of taxa. Philippe el al (1994) 
have showed that by eliminating rapidly evolving species 
from the analysis, discrepancies between molecular and 
traditional phylogeny partly disappear and bootstrap val¬ 
ues rise at some nodes. Since we only intended to show 
the placement of the dicyemids firmly within the triplo¬ 
blasts, we chose taxa representing a broad spectrum of 
eukaryote phyla without regard to consistency with the 
traditional view of triploblast phylogeny. 

There have been opposing views on the role of the di- 
cyemid mesozoa in the story of metazoan evolution. Some 
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Figure 3. The consensus tree obtained using the maximum-parsimony algorithm with bootstrap resam¬ 
pling {DNAPARS, SEQBOOT. and CONSENSE programs of PHYLIP package, version 3.5c), showing the 
phylogenetic position of the dicyemids among 28 representative eukaryolic taxa. The percentage of 500 
parsimony replicates is shown at the node the value is supporting. The tree was rooted by using Saccharomyces 
cerevisiae as an outgroup. Branch lengths are proportional to the scale given in number of substitutions (a 
total of 2807). T his tree is different from the three most parsimonious trees in the positions of Crassostrea 
gigas and Mohniformts molinifnmus within the assemblage of the coelomate triploblasts. 


authors have proposed that the dicyemids arc a missing 
link between unicellular organisms and multicellular an¬ 
imals (Dodson, 1956; Hyman, 1959; Lapan and Morow- 
itz, 1974; Ohama el al ., 1984), while others have claimed 
that they are an animal group degenerated as a result of 
parasitism (Nouvel, 1948; McConnaughcy, 1951; Stun- 
kard, 1954; Ginetsinskaya, 1988). The phylogenetic trees 
inferred from comparisons of nucleotide sequences of 5S 
rRNA suggested that the dicyemids emerged first among 
the metazoa examined and that triclad flatworms, nem¬ 
atodes, cnidarians, and sponges followed, in that order 
(Ohama el al., 1984, Hori and Osawa, 1987). This sug¬ 
gestion does not, however, accord with the present result 
and the previous inferences about metazoan phylogeny 
based upon I8S and 28S rDNA sequences (Field el al., 
1988; Christen el al., 1991; Wainright el al., 1993; Ko- 
bayashi el al., 1993). Discrepancies are partly ascribable 
to differences in the methods used to infer phylogenetic 
relationships. In contrast to the 18S and 28S rDNA trees 
reconstructed by the neighbor-joining, maximum-parsi¬ 
mony, and maximum-likelihood methods, the above 5S 


rRNA trees have been reconstructed by unweighted and 
weighted pair group methods using arithmetic averages 
(UPGMA and WPGMA, respectively), which are valid 
under the assumption that rates of nucleotide substitution 
are constant among taxa analyzed (Sokal and Mitchener, 
1958). However, the essential point is that the 5S rRNA 
is too small to contain signal sufficient to allow precise 
inference of phylogenetic relationships. Because of large 
standard errors, sequential orders of branching of the di¬ 
cyemids, flatworms, nematodes, cnidarians, and sponges 
shown i n the above 5S rRNA trees appear to be statistically 
insignificant. Recently Halanych (1991) analyzed the se¬ 
quence data of 5S rRNA with the maximum-parsimony 
method. The phylogenetic tree obtained was inconsistent 
with phylogenies based on 18S and 28S rDNA data, and 
few nodes in the tree were supported by bootstrap value 
at a significant level. 

The present results do not appear to support the prop¬ 
osition that the dicyemids are a truly primitive group 
linking unicellular organisms with multicellular metazoa. 
Instead, our results favor the view that the dicyemids are 
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Figure 4. Maximum-likelihood tree showing the phylogenetic position of the dicyemids among 28 rep¬ 
resentative eukaryotic groups. The tree was obtained using the fastDNAML algorithm with a transition/ 
transversion ratio of 1.48, which gave the best maximum-likelihood score (In likelihood = -11909.45548). 
Branch lengths are proportional to the scale given in substitutions per sequence position. Because of the 
long computation time (more than 12 h per replication), bootstrapping was not performed. 


degeneratively simple animals descended from a more 
complex triploblastic ancestor. Recent close observations 
of dicyemid development (Furuya el al. t 1992b, 1994) do 
not contradict the present inference; spiral cleavage, a de¬ 
terminative mode of cell division, and the formation of 
stereoblastula-like structure through epiboly, as seen in 
the development of vermiform and infusoriform embry os, 
are reminiscent of flatworms. Mvxozoans arc also re¬ 
garded as an extreme example of the degeneration char¬ 
acteristic of parasitic evolution (Smother el al., 1994). At 
present we can say little about the ancestor from which 
the dicyemids were derived. Molecular analyses that in¬ 
clude other lower turbellarian groups (for example, the 
Catenulida and the Nemertodermatida) and the Ortho- 
nectida, a group tentatively included in the Mesozoa, will 
provide further information for understanding the phy¬ 
logenetic position of the dicyemids. 
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